arXiv:1509.06307vl [nucl-th] 21 Sep 2015 


Toward a simultaneous description of Raa and V2 for 
heavy quarks 

S. K. Das^’^, F. Scardina S. Plumari^’^ and V. Greco^’^ 

^ Department of Physics and Astronomy, University of Catania, Via S. Sofia 64, I- 95125 
Catania, Italy 

^ Laboratori Nazionali del Sud, INFN-LNS, Via S. Sofia 62, 1-95123 Catania, Italy 
E-mail: Scintosh@lns.infn.it 


Abstract. The two key observables related to heavy quarks that have been measured in RHIC 
and LHC energies are the nuclear suppression factor Raa and the elliptic flow V2- Simultaneous 
description of these two observables is a top challenge for all the existing models. We have 
highlighted how a consistent combination of four ingredients i.e the temperature dependence of 
the energy loss, full solution of the Boltzmann collision integral for the momentum evolution 
of heavy quark, hadronization by coalescence, then the hadronic rescattering, responsible to 
address a large part of such a puzzle. We have considered four different models to evaluate the 
temperature dependence of drag coefficients of the heavy quark in QGP. All these four different 
models are set to reproduce the same Raa as of the experiments. We have shown that for the 
same Raa, the V2 could be quite different depending on the interaction dynamics as well as 
other ingredients. 


1. Introduction 

Theoretical calculations predict that hadronic matter at high temperatures and densities 
dissolves into a deconfined state of quarks and gluons - called Quark Gluon Plasma (QGP).The 
experimental efforts at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) 
energies is aimed at to create and characterized the properties of QGP. The heavy hadrons 
(hadrons which contain at least one heavy quark, mainly c and b ) constitute a unique probe of 
the QGP properties, because they are produced in the early stage of the collisions and they are 
the witness to the entire space-time evolution of the system. 

The two key observables related to heavy quarks that have been measured in RHIC and LHC 
are the nuclear suppression factor Raa which is the ratio between the the pT spectra of 

heavy flavored hadrons (D and B) produced in Au-|-Au collisions with respect to those produced 
in p+p collisions, and the elliptic flow V 2 mil], which is a measure of the anisotropy in the 
angular distribution. Several theoretical efforts have been made to study the Raa and the V2 
measured in experiments within the Fokker Planck (FP) approach O [6l [3 [HI [9l [TO] and the 
relativistic Boltzmann approach imiiziiiaEiiniiiBi- However all the approaches shown some 
difficulties to describe both the Raa and V2 simultaneously. 

2. Results 

To address the simultaneous description of Raa and V2, we start with the time evolution of 
Raa and V2 to know how they develop during the expansion of the QGP. To study the time 


evolution of Raa and ^ 2 , we have solved the Fokker Planck (FP) equation stochastically in terms 
of the Langevin equation, for detail we refer to our earlier work [TO]. The elastic collisions of 
heavy quarks with the bulk has been considered within the framework of pQCD. The divergence 
associated with the t-channel diagrams due to massless intermediate particle exchange has been 
regularized introducing the Debye screening mass mn = dvra^T with the running coupling 
For the bulk evolution we are using a transport bulk which can reproduce some of the gross 
features of the bulk naiis]. At RHIC energy, Au + Au at ^/s = 200, we simulate the fireball 
with initial temperature in the center is Tj = 340 MeV and the initial time is r* = 0.6 fm/c. 
for the detail of the initialization, we refer to our early work. Initially the charm quark are 
distributed according to the charm production in pp collisions |2D] . 




Figure 1. Time evolution of Raa- Figure 2. Time evolution of V 2 . 

As shown Figlll the Raa is such which is develop at the very early stage of the evolution and 
get saturated within 3-4 fm where as the V 2 (in [2|) is develop at the latter stage of the evolution. 
Initially the bulk have zero v2. First the bulk will develop its own V 2 and then it will transfer it 
to the heave flavor. Hence the Raa is sensitive to the magnitude of the drag coefficient at the 
early stage of the evolution i.e at R where as the V 2 is very sensitive to the magnitude of the 
drag coefficient at the latter stage of the evolution i.e at Tc . This highlights the T-dependence 
of the drag coefficient may play a significance role for a simultaneous description of Raa and 
V 2 as they are sensitive to different range of T. To study the impact of T-dependence of drag 
coefficients on heavy quark observable, we have consider four different model to calculate the 
drag and diffusion coefficients of heavy quark in QGP. 

Model-I (pQCD): The elastic collisions of heavy quarks with the bulk has been considered 
within the framework of pQCD. 

Model-II (AdS/CFT): In this case we have considered the drag force from the gauge/string 
duality i.e. AdS/CFT [21] which reads Tconf = C ; where C = 2.1 ± 0.5.(also see [22] i 
Model-Ill: In the third case, we have employed a quasi particle model (QPM) model [23l [Ml IM] 
with T-dependent quasi-particle masses, niq = l/3g‘^T‘^, rUg = 3/4g'^r^, plus a T-dependence 
background field known as bag constant, tuned to the thermodynamics of the lattice QCD. Such 
a fit lead to the coupling, g‘^{T) = T/Ts=0.57. 

Model-IV {aQPM{T),mg = rUg = 0) In this case, we are considering a model where the light 
quarks and gluons are massless but the coupling is taken from the QPM model discussed in 
Model III. 

This fourth case has been considered to have a drag which decreasing with T as obtained 
in the T-matrix approach [61126] .The T-dependence of the drag coefficients have been shown in 
Fig [3] obtained within the four model discussed above at pp = 100 MeV. These are the rescaled 






drag coefficients which can reproduced the same Raa as of the experiment at RHIC energy. 





Figure 3. Drag coefficients 
as a function of temperature. 


Figure 4. Raa as a function Figure 5. V 2 as a function of 
of pt- Pt- 


In Fig 0] we have shown the Raa as a function of pT for the four different models obtained 
within the Langevin dynamics at RHIC energy. The V 2 for the same Raa has been plotted in 
Fig [5] for all models as a function of pT- We found that for the same Raa, the V 2 build up can 
be quite different depending on the T-dependence of the interaction. Larger the interaction at 
Tc larger is the V 2 [lU]. Similar effect has also been found in light quark sector [271 [2H] . This 
suggests, the correct temperature dependence of drag coefficient has a significance role for a 
simultaneous reproduction of Raa and V 2 - 



Figure 6. Raa as a function 
of Pt at RHIC energy for QGP 
and QGP+hadronic phase. 


Figure 7. V 2 as a function of 
Pt at RHIC energy for QGP 
and QGP+hadronic phase. 


Figure 8. Raa vs V2 at pt 
1.5 GeV. 


Recently it has been shown m that in case of the full solution of the Boltzmann (BM) 
integral i.e. without the assumption of small collisional exchanged momenta, leads in general 
to a large V 2 than that of the Fokker Plank case which is a approximation of the Boltzmann 
equation under the assumption of small exchanged momenta. Also heavy quark hadronization 
by coalescence further enhancement the V 2 [61129]. 

To study the role of the hadronic phase (HP) on the heavy quark observables, we study 
the propagation of D mesons in the hadronic medium consist of pions, kaons and eta. The 
elastic interaction between the D meson with the hadronic matter has been treated within the 
framework of chiral perturbation theory [30] . In Fig[6]we have shown the impact of the hadronic 
medium on Raa which is almost unnoticeable. This is because the Raa develop at the early 
stage of the evolution and get saturated due to radial flow within 3-4 fm. The impact of the 
hadronic medium on V 2 has been shown in Fig[7| The hadronic medium further enhance the V 2 
around 20%. In Fig [8] we have shown how the V 2 build up depending on the T-dependence of 
the drag coefficients for the same Raa- Then the V 2 get a boost from the Boltzmann equation 















(BM) in terms of evolution over the Fokker Planck equation. The V 2 further enhanced due 
to hadronization by coalescence. Finally the V2 get the hadronic phase boost where the Raa 
remain the same. 

3. Summary and outlook 

In summary, we have shown how the V2 build up for the same Raa depending on the T- 
dependence of the drag coefficients which is key for a simultaneous description of Raa and V2- 
We have also highlighted how the V 2 get a boost from the Boltzmann dynamics to study the 
momentum evolution of heavy quark and from hadronization by coalescence. Then we have 
shown the effect of the hadronic medium on Raa and V2- The impact of radiative processes 
on heavy quark observables will also be investigated in an upcoming study along with all these 
four ingredients within a single framework . It will also be interesting to study the role of 
pre-equilibrium phase m on Raa and V2 relation. 
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